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Abstract. Motivated by G. H. Hardy's 1939 results [4] on functions orthog- 
onal with respect to their real zeros A n , 71 = 1, 2, ... , we will consider, within 
the same general conditions imposed by Hardy, functions satisfying an orthog- 
onality with respect to their zeros with Jacobi weights on the interval (0,1), 
that is, the functions f(z) = z"F(z), v SE R, where F is entire and 
r 1 

f(X n t)f(\ m t)t a (l - tfdt = 0, a > -1 - 2v, P > -1, 

^0 

when n 7^ m. Considering all possible functions on this class we are lead to 
the discovery of a new family of generalized Bessel functions including Bessel 
and Hyperbessel functions as special cases. 



■hi 



1. Introduction 

In his 1939 paper [4], G. H. Hardy proved that, under certain conditions, the 
only functions satisfying 

»i 

f(\ n t)f{\ m t)dt = 





if m 7^ n, are the Bessel functions. 

With a view to extend Hardy's results to the g-case, it was observed in pQ that 
a substantial part of Hardy's argument could be carried out without virtually any 
change, when the measure dt is replaced by an arbitrary positive measure du(t) on 
the interval (0, 1) and functions satisfying 

(1.1) / f(X n t)f(X. m t)du(t) = 







if m n, are considered. 
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In particular, the completeness of the set {/(A ra i)} holds for such a general 
orthogonal system. Furthermore, there exists an associated Lagrange-type sampling 
theorem for functions defined by integral transforms whose kernel is defined by 
means of / pQ . 

It is however not possible to identify, under such a degree of generalization, the 
whole class of functions satisfying (jl.ip . Therefore, we raise the question: Given a 
specific positive real measure on the interval (0, 1), what are the resulting orthogonal 
functions? In the case where this measure is the d q x arising from Jackson's q- 
integral, it was shown in pQ that the corresponding functions are the Jackson q- 
Bessel functions of the third type. In this note we will answer this question in 
the case where du(t) = t a (l — tfdt is the measure associated with the orthogonal 
Jacobi polynomials on the interval (0, 1). That is, we will characterize in some sense 
the functions f(z) — z v F(z), where v 6 K, F is entire, which satisfy the equality 

(1.2) f f(\ n t)f(\ m t)t a (l - tfdt = 0, a > -1 - 2v, (3 > -1, 
Jo 

with membership of the classes A or B being specified in the following definition 
(the same general restrictions are imposed in [4] and [T]). 

Definition 1. Let v,a el be such that 2v + a > -1 and f{z) — z"F(z). The 
class A is constituted by entire functions F of order less than two or of order two 
and minimal type having real positive zeros X n > 0, n = 1,2, ... with ^ A~ 2 < 
oo situated symmetrically about the origin. The class B, in turn, contains entire 
functions F(z) which have real but not necessarily positive zeros X n , n = 1,2, ... 
with ^2 X~ < oo and of order less than one or of order one and minimal type with 
F(0) = 1. 

In this paper we will deal mainly with functions of the class B. Namely, we will 
show that these functions satisfy the Abel-type integral equation of the second kind 
[7] and will study their important properties and series representations. However, 
in the last section we will derive similar integral equation for the class A and will 
also give its explicit solutions in terms of the series. 

We will appeal in the sequel to the theory of the Mellin transform [BJ, ]T\, [§]• 
As it is known the Mellin direct and inverse transforms are defined by the formulas 

/>OC 

(1.3) f*(s) = / f(x)x s - 1 dx, 

Jo 

i p'y+ica 

(1.4) f(x) = — \ r(s)x- s ds, s = 7 + xt, x > 0, 

where integrals (1.3)- (1.4) exist as Lebesgue ones if we assume the conditions 
/ G Li(R + ; x~ f ~ 1 dx), f* € ii((7~ ioo, 7-M00); dr), respectively. However functions 
from the classes A or B can have a non- integrable singularity at infinity and the 
classical Mellin transform (1.3) generally does not exist. Therefore owing to [8], Ch. 
4 we introduce the Mellin transform for distributions from the dual space MX ^ 
to Zemanian's testing -function space M 7i 72 into the space of analytic functions 
in the open vertical strip f2y := {s S C : 7 X < Res < 72} by the formula 



(1.5) 



f*(s) := {f (x^x*- 1 ), seQf. 
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In this case the inversion integral (1.4) is convergent in D'(R+), i.e. for any smooth 
function 9 6 2?(R+) C M 7 7 with compact support on R+ it holds 



(1.6) (f(x),e(x)}= hm(— / f*(s)x- s ds, 9{x)), 

r—xx) ZTTl 



■y—tr 



2. The Hardy type integral equation 

Suppose that f £ B and satisfies (|1.2j) . Let 



l 

[/(A„i)] 2 i Q (l - tfdt, a>-l-2v, > -1. 

If a„(z) stands for the Fourier coefficients of the expansion of f(zt) in terms of the 

_ i 

orthonormal basis {A n 2 f{X n t)}, then it is possible to show, following the arguments 
in [4], that 

(2.1) «»(*) - ^ / /(-)/(Wrd - - /.(A^'-A,,) - 

In the sequel we will essentially continue to go along the lines of [4], but since 
considerable number of adaptations have to be made in order to deal with the 
presence of the Jacobi weight, we find it better to provide the proof. 
Using the Parseval identity for inner products we get 



(2.2) 



-1 oo 

/ f(zt)f((t)t a (l - tfdt = V a n (z)a n ((). 
Jo 



Theorem 1. Let f e B and satisfy with a > -1 - 2v, (3 > -1. 

Then it satisfies the integral equation 

a I u v+a+1 (z-uff(u)du = (az + l) [ u v+a {z - uf f{u)du 
Jo Jo 

(2.3) +z v+a+l3+1 f{z)A, 

where a = F'(0) and A = —B(2v + a + 1, /3 + 1), B denotes the Beta- function [3]. 
PROOF. Substituting (2.1) into (2.2), we obtain, after some simplifications, 

(2.4) f 1 f(zt)f((t)t a (l - tfdt = -f(z)f(O q{z) ~ qi0 



where 



= 2. 7777 



1 



z — A„ Ar, 



- {/'(A„)} 2 
Letting £ — > in (2.4), the result is 

t t» +a f(zt)(l - tfdt = - /(z)g(2) . 
Jo z 

The change of variables u — zt gives 

(2.5) / U v+a (z-uff(u)du = -z v+a+l3 f(z)q{z). 
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Now, when z is small enough, f(z) ^ z v and q(z) ~ zq'(0). Therefore, when z — » 
we have 



2i/+a 



(z - u^du - -z 2t/ + Q+/3 +Y(0) 



and, as a consequence, 

q'(0) = -B(2u + a + l,/3+l). 
Rewriting (2.4) in the form 

t 2u+a F(zt)F((t)(l - tfdt = -F(z)F(Q q ^~ q ^ 

and differentiating in £ we obtain 

t 2v+a+1 F'(Ct)F(zt)(l-tfdt 



Setting £ = 0, it gives 

F'(0) / f +a+1 (l - tff(zt)dt = [q(z)(F'(0)z + 1) - A?] . 



Jo 

So, if a = F'(0), the change of variable zt = u yields 



a l u 

o 



u+a+l 



(z - uff{u)du = -z u+a+l3 f{z) [(za + l)q(z) - Az] 



Combining with (2.5) we get the integral equation (2.3). □ 

3. The Hardy-type theorem 

In this section we will prove the Hardy type theorem [4] , which will characterize 
all possible orthogonal functions with respect their zeros from the class B under 
the Jacobi weight. Indeed, we have 

Theorem 2. Let f satisfy conditions of Theorem 1. Then 

oo 

(3.1) f(z) — const. z v a n z n , 
where 

(3 . 2 ) a ,, = M±llfn M; 

T{fi + n) (/x + f3 + l)j - (/x)j 

-ffere /i = 2z/ + a + 1, T(z) is Euler's Gamma-function, (b)j is the Pochhammer 
symbol [3] and the empty product is equal to 1. Moreover, series (3.1) represents 
an entire function of the order p = < 1, when (3 > — 1. In particular, the case 
f3 = drives at the Hardy solutions in terms of the Bess el functions (cf. [4], p. 43/ 

Proof. Making again an elementary substitution u = zt and cutting the mul- 
tiplier z u + a +P+ 1 we write (2.3) in the form 

az [ t v+a+1 (l-tff(zt)dt={az+l)[ t u+a {l - tf f(zt)dt 
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(3.3) -B{2v + a + l,0 + l)f{z). 

Hence we observe that (3.3) is a second kind integral equation containing two 
Erdelyi-Kober fractional integration operators with linear coefficients (see [7] , Ch. 
3). Seeking a possible solution in terms of the series z v ^2^ =Q a n z n where a n ^ 
0, eto = 1 and no two consecutive a n vanish (see [4], p. 43) we substitute it in 
(3.3). After changing the order of integration and summation via the uniform 
convergence, calculation the inner Beta-integrals, and elementary substitutions we 
come out with the equality 

oo 

a a n -iz n+y [B(n + n, + 1) - B(p + n- 1,0 + 1)] 



(3.4) 



OO 

n=l 



[B(ji + n,P + l)- B(n, + 1)}, n = 2v + a + l. 



Hence equating coefficients of the series and taking into account that B(p + n, + 
1)-B{fx,0 
relations 



l)/0,neN,/i>0, > — 1 we obtain the following recurrence 
B(fj, + n,0 + l)-B(fj, + n- 1,0 + 1) 



(3.5) 



B(n + n,p + l)-B(p,p+l) 
B(» + j,0+l)-B(fM + j -1,0+1) 



B(fx + j,0 + l)-B(^,0+l) 

where n e No and the empty product is equal to 1. Moreover, invoking the defini- 
tion of the Pochhammer symbol (b)k = ^n?) i the reduction formula for Gamma- 



_ r(fc)r(c) 



we 



function T(z + 1) = zT(z) and representing Beta-function via B(b, c) — ^(b+c) 
write (3.5) in the form (3.2). 

Further, it is easily seen from (3.5) and asymptotic formula for the ratio of 
Gamma-functions [3] that 

B(p + n,0+l)- B(p +n- 1,0 + 1) 



B(p + n,0+l)- B(p,0+1) 
n-1) \a>\{0+!) 



T(p + n + 0+l) B(p, 0+1)- B(ji + n,0+ I) 

= O (n~ (/3+2) ) -> 0, n -> oo, > -1. 

Thus F(z) — X)ri°=o a « z ™ i s an entire function. Let us compute its order p. To do 
this we appeal to the familiar in calculus Stolz's theorem. Hence we find 



p = limsup„_ 



= lim 



nlogn 

logKl- 1 

log 71 



= lim 



(n + 1) log(n +1) —n log n 



log | a. 



'i+i 



log|a„ 



f 



< 1, 0> -I. 



logn^+ 2 + 2 

Thus same arguments as in [4] guarantee that all our solutions belong to the class 
B. In particular, for = we easily get from (3.2) that 



n 



(apY 



{ji + n) (n + l)j - (jj,)j T(p + n)n\ 



n = 0,l,. 
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Then the series representation for the Bessel function [3] yields 
f(z) = const. z~ a/2 J p _i (c- 2 ^ 2 ) ! c 2 = -4a//, 
which slightly generalizes Hardy's solutions [4], p. 43. Theorem 2 is proved. □ 

4. Properties of solutions and their particular cases 

First in this section we will apply the Mellin transform (1.5) to the integral 
equation (2.3), reducing it to a certain functional equation. Then we will solve this 
equation by methods of the calculus of finite differences [5] to obtain the value of 
the Mellin transform for solutions (3.1). 

In fact, returning to (3.3) we consider z = x e M + and we apply through the 
Mellin transform (1.5) taking into account the operational formula 

(/(zt),^- 1 )^" 8 /*^). 

Then using values of the elementary Beta-integrals we come out with the following 
homogeneous functional equation 

aB{y + a- s + l,/3+ l)/*(s + 1) = aB{y + a-s,(3+ l)f*(s + 1) 

(4.1) +B{v + a-s + l,f3+ l)/*(s) - B(2u + a +l,p+ l)/*(a), 

where s is a parameter of the Mellin transform (1.5) such that s + l£fi/. Hence 
denoting by h(s) = B(v + a — s + 1,/3+1) and invoking the condition o ^ (see 
@], p. 43) we rewrite (3.2) as 

(4.2) t(s+l) = H(s)r(s), 
where the kernel H(s) is given by 

1 h(s) - h{-v) 



(4-3) H(s) 

a h(s) — h{s + 1) 

We can simplify H(s) appealing to the basic properties for Beta-functions |3j. After 
straightforward calculations we get finally 

(44) H{8)= 1 . (s-u-ams)-h(-u)) 

a(p + l) h(s) 

Meanwhile, 

" (-/3)„ 



h(s) - h(-v) = / (1 - tf{t u+a - s - t 2v+a )dt = V 

JO n 



! 

n=0 



1 1 



where 



X(s) = E 



1 2v + a + n-\ 



n\(2v + a + n + l)(v + a + n + 1 - s) ' 

n— ' 

It is easily seen that x(s) and h(s) have the same simple poles at the points s = 
i/ + a+n + l, n = 0, 1, . . . .As a consequence, taking into account the expression of 
h(s) in terms of Gamma-functions, we find that H(s) is meromorphic with simple 
poles s = is + a + (3 + n + 2, ne No- Furthermore, via the asymptotic behavior of 
the Beta-function [3], we get H(s) = 0(s fJ+2 ), s — > oo. 
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The functional equation (4.2) can be formally solved by methods proposed in 
[5], Ch. 11. So it has a general solution 

(4.5) /*(*)= u(s)exp^logff(z)Az^ , 

where u>(s) is an arbitrary periodic function of s of period 1 , J2l means the operation 
of general summation [5) with a given constant c and the function log_ff(s) is 
summablc in this sense (cf. [5] Ch. 8). Clearly that it may be necessary to make 
suitable cuts in the s plane in view of the possible multi- valued nature of the right- 
hand side of (4.5). 

Let us consider interesting particular cases of the solutions (3.1) when (3 = 
fc £ No. Indeed, putting (3 = k in (3.2) and using the definition of Pochhammer's 
symbol we obtain 

(n + k + i)j = r(/i + fc + j + i) (mWj+i _., 



_ + + j + 1) ...(// + fc + j) 1 _ jPfc(j) 

(M)fe+i Mfc+i ' 

where we denote by Pk(j) = (j — Oii) ■ ■ ■ (J — ctk), a t £ C, i = 1, . . . , k, a polynomial 
of degree fc with respect to j and evidently Pk(j) > 0, j G N. Substituting the 
above expression into (3.2) it becomes 

(Art ( a (fc + l)fc) fc+1 )" An 1 , cN 

(4 - 6) a - = r( M +n)n! nnj^' n ' e ° 

and the empty products are equal to 1. Consequently, solutions (3.1) represent the 
so-called hyperbessel functions [7], Ch. 19 of the order k + 1, namely 

/(z) = const. z y 

(4.7) 

Xo-Ffc+i (2z/ + a + 1, 1 -ai,..., l-a fc ; a(fc + l)(2i/ + a + l) k+1 z) , fc € N . 
As it is known [7j functions 

F(z) = F fc+ i (2i/ + a + 1, 1 - ai, . . . , 1 - a k ; a(k + l)(2v + a + l) k+1 z) 
in (4.7) satisfy the following k + 2-th order linear differential equation 

(4.8) ^ (z^- +2v + a\]J (z^- - a^j F - a(k + l)(2u + a + l) k+1 F = 0. 

Remark 1. We point out that in [2] the authors considered other generaliza- 
tions of the Bessel functions which also satisfy higher order differential equations. 

Remark 2. In particular, for k = we easily get Hardy's solutions in terms 
of the qFi - functions and strictly prove their orthogonality (1.2) by using the 
corresponding second order differential equation (4.8). However, even for k = 1 the 
direct proof of the orthogonality property (1.2) for solutions (4.7) is a difficult task 
since we deal in this case with the third order differential equation. As we aware, 
the orthogonality property for the hyperbessel functions with respect to their zeros 
is yet unknown. So, following conclusions of Theorem 2 we conjecture here this fact 
as well as the orthogonality for all solutions (3.1) of the class B. 
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5. Orthogonal functions of the class A 

In the case / G A we recall again the arguments in [4] to write accordingly, the 
formula for Fourier coefficients a n (z) of f(zt) (see (2.1)) as 

(5.1) a„(z) 



f\\ n ) z 2 -A 2 " 
We have 

Theorem 3. If f e A and satisfy l\1.2\i with a > -1 - 2u, u € M, /3 > -1, 
£/ien £/ie integral equation holds 

a f u u+a+2 (z-uff{u)du 



(5.2) = {az 2 + 2) u u+a {z-uff(u)du + z v+a+ ^ +1 f(z)A, 

Jo 

where a = F"(0) and A = -2B(2v + a + 1,(3+1). 

Proof. As in the proof of Theorem 1 we substitute (5.1) in (2.2) and we get 
the equality 

(5.3) f(zt)f(tt)ni - tfdt = ~f(z)f(o q{z z lip\ 



where now 



* A \ 2 



2 \ 2 \ 2 



Letting £ — > in (5.3), we obtain 

/' f /(zt)t a (l - tfdt = JWf- z \ 
Jo z 

which yields after the change u = zt 

(5.4) / u" +a {z-u) l3 f{u)du^-z u+a+l3 - 1 f{z)q{z). 



When z is small, f(z) ^ z v and q(z) ~ ^-g"(0). Therefore, as z — > we get 

2 f u 2 "+ Q (z - ufdu ~ -2 2l/ + a+/3+1 g"(0). 
Jo 

Thus q"(0) = ~2B(2v + a + 1, (3 + 1) = A. If we rewrite (5.3) in the form 



2 z-2 



^ +Q F(zt)F(a)(l - tfdt = -F(z)F(C)—^ — 
then after differentiation with respect to £ we get 
t 2v+a+1 F'(tt)F(zt)(l-tfdt 
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F"(0) [ t v+a+2 {\ - tff(zt)dt = -f(z) 
Jo 



Differentiating again in £ and setting £ — 0, we deduce after some simplifications 

< F"{Q)q{z) A | 2q{zy 
z 2 z 2 z 4 

Letting a = F" (0) and making the change zt = u we obtain the equation 

a f u v+a+2 {z - uY 3 f{u)du = -z u+a+fj -\f{z) [(z 2 a + 2)q(z) - Az 2 
Jo 

Taking into account (5.4) we finally obtain the integral equation (5.2). 
An analog of the Hardy theorem for the class A is 
Theorem 4. Let f satisfy conditions of Theorem 3. Then 



□ 



(5.5) 

where 



f(z) — const. z u a2nZ 2n , 



71=0 



a([3 + l) 



n (/3+2( ^ 2 -!: 1))( ^-^ r , = o.i.. . 



(At + P + l)2j - (ji)2j 
ji = 2v + a + 1 and the empty product is equal to 1 . 

PROOF. The substitution u = zt and the value of A reduce equation (5.2) to 

^2 



az 



f 1 t v+a (t 2 -tff(zt)dt 
Jo 



(5.6) 



= 2 



f 1 t v+a {\ - tff(zt)dt - f{z)B{2v + a + 1, (3 + 1) 
Jo 



Since the case / £ A presumes the following series representation 



f(z) = const. z v ^ a,2„; 



In 



n=0 



with a2n ^ for any n <G No, we substitute it into (5.6), change the order of inte- 
gration and summation, and calculate the inner Beta-integrals. Therefore denoting 
by [i = 2v + a + 1 we obtain 

oo 

a a 2{n -i)Z 2n [B(n + 2n, (3 + 1) - B(p + 2(n - 1), [3 + 1)] 



n=l 



2^> 2 „z 2n [i?( M + 2n,/?+l) - /?+!)]. 



n=l 



Hence equating coefficients of the series the following recurrence relations appear 

_ a B(n + 2n, + 1) - gQx + 2(n - 1), g + 1) 

^-^(n-i) B(/i + 2 „,/J+l) - B(jm,0 + 1) 

'a\ n A S 0" + 2 i> /9 + 1) - + 2(7 -!),/?+ 1) 



(5.7) 



(1) n 



B(/x + 2i,/3 + l)- J B(/x,/3 + l) 
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where n = 0, 1, 2, ... , and the empty product is equal to 1. In the same way as 
(3.2) it can be written in the form 



l2n 



ajfi + 1) \ n fr (/3 + 2Qx + 2j - 1))(m) 2(j --i) 

11 ( M + ( 3 + l) 2 ,-( M )2j ' ^ U ' ' 



□ 

Putting j3 = in (5.7) by straightforward calculations functions (5.5) become 

f(z) — const. z^~Jh_ 1 ^(-a^) 1 / 2 ^^ , 

which slightly generalizes Hardy's solutions [4], p. 42. 

It is easily seen that the series in (5.5) is an entire function. Finally we compute 
the order of these solutions. By similar arguments as in Theorem 2 we have 

2nlog2n , 2(n + 1) log2(n + 1) - 2nlog2n 

p = hmsup^^- : pj- = hm — . — : — — 

log|a 2n | ™^°° log|a 2 („+i)| 1 -log|o2n| 

n^oo logn^ 2 (3 + 2 
Therefore all solutions belong to the class A. 

Remark 3. Our final conjecture is that solutions / S A of the Hardy-type 
integral equation (5.2) are orthogonal with respect to their own zeros for the Jacobi 
weight. 
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